Inertial and fluctuational effects on the motion of a Bose superfluid vortex 
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We study the motion of a vortex under the influence of a harmonic force in an approximately two 
dimensional trapped Bose-condensed gas. The Hall-Vinen-Iordanskii equations, modified to include 
a fluctuational force and an inertial mass term, are solved for the vortex motion. The mass of the 
vortex has a strong influence on the time it takes the vortex to escape the trap. Since the vortex 
mass also depends on the trap size we have an additional dependence on the trap size in the escape 
time which we compare to the massless case. 
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INTRODUCTION 

Quantum vortices appear in a large range of phys- 
ical systems ranging from laboratory condensed mat- 
ter systems like superconductors [ij and neutral super- 
fluids (^He superfluid superfluid ^He, and cold atomic 
gases [3|), to exotic excitations playing a role in the early 
universe y. Understanding the motion of these vortices 
is key to understanding the properties of such systems. 

The motion of superfluid vortices has been studied 
since the 195O's0. Their dynamics has often been de- 
scribed by some form of Hall-Vinen-Iordanskii (HVI) 
equation^, which essentially describes the different 
forces acting on a vortex. However until recently the form 
and even existence of both the forces 0, Q and the vor- 
tex mass [9| ha ve been controversial. Recently Thompson 
and Stamp [lO| derived an equation of motion for a vortex 
in a boson superfluid, starting from a low energy effec- 
tive field theory for the superfluid. In the limit of low 
frequency motion, which turns out to be a semiclassical 
limit, this equation of motion reduced to a modified HVI 
equation, in which both an inertial term and a Markovian 
noise term are added to to the original HVI equation. 
However it was also found [lo[ that at higher frequen- 
cies, strong departures occurred from the HVI equation, 
and the vortex coordinate then obeys a different equa- 
tion which shows strong retardation effects, with highly 
non- Markovian correlations. 

In this paper we show how the HVI equations, modi- 
fied to include the inertial and noise terms, can be used 
to discuss the dynamics of a single vortex in a cold BEC 
gas. At the temperatures and frequencies so far em- 
ployed in experiments on such gases, it turns out we are 
in the well within the semiclassical regime. Vortices in 
trapped Bose-Einstein condensates were first created in 
the lab by Matthews et al. [Tl[ , and since then single vor- 
tices and lattices of vortices in BECs have been studied 
in various contexts 12 - 1^ . The motion of a single vor- 
tex in a trapped BEC is the subject of this paper. As 
predicted by Rokhsar [lij and observed later observed 
experimentallvflll [l^ single vortices created in an opti- 
cal trap spiral out from the centre of the trap until they 



decay at the edge of the condensate. 

There have been fairly extensive theoretical studies of 
vortex dynamics in BECs; however these have not in- 
cluded an inertial term. Jackson et al. 17| and Svidzinsky 
and Fetter [l^ used the Gross-Pitaevskii equation to cal- 
culate the precession velocity of the vortex. Fedichcv and 
Shlyaponikov [l^ used a version of the HVI equations 
with a term describing dissipation, derived by consider- 
ing scattering of the single particle excitations from a sta- 
tionary vortex core and a backgroud superflow due to the 
interaction of the vortex with the boundary of the con- 
desate. More recently Jackson et al.^^ performed nu- 
merical simulations, and interpreted the results in terms 
of an HVI equation with no inertial mass term or back- 
fiow. Duine et al. [2lj used a stochastic form of the 
Gross-Pitaevskii equation at finite temperature to derive 
a purely dissipative equation for the motion of the vortex, 
including a thermal noise term. 

In what follows we study the motion of a vortex in a 
trap using the form of the two-dimensional HVI equation 
derived by Thompson and Stamp [lo| in the semiclassical 
limit. We assume the vortex is in a harmonic potential 
caused by the trap. We solve these equations of motion 
for the motion of the vortex and calculate the time it 
takes the vortex to escape the trap. We discuss how the 
escape time depends on the radius of the trap and on 
the temperature, and we discuss the role that the vortex 
inertial mass plays in these results. We also compare the 
results to previous work. 



EQUATION OF MOTION 

A quantum vortex can be described either by look- 
ing at the A^-particle wave-function of the system, or 
by defining a vortex reduced density matrix p(r,r';t) = 
(r|/5(t)|r'), in which the other degrees of freedom of the 
superfluid have been averaged over. If we wish to find 
an equation of motion for the vortex itself, we are then 
obliged to derive this from the dynamics of p{r,r';t). 
This density matrix is conveniently rewritten in terms 
of a 'centre of mass' variable R = (r + r')/2 and a quan- 
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turn fluctuation variable ^ = (r — r')/2. In the 'classical 
limit' where the characteristic frequency f2 of the vor- 
tex dynamics is low (such that Ml <^ kT, where T is 
the temperature), the quantum fluctuations (^it) become 
negligible, and we then expect [l^l a set of modified HVI 
equations to be valid for what is now a semiclassical vor- 
tex coordinate R(t). In a 2-dimensional Bose superfluid 
these take the form 

. , d'^R. „ dR, „ dR „ , . „ / N N 
M„^-pKzx— +i?o(r)— =Ffluc(i) + Fap(i) (1) 

in a frame of reference where v„ = Vs = (ic, botht 
he normal current and supercurrents are zero; for the 
more general case where they are arbitrary, see ref . (loj ) . 
Here Fap is the applied force on the vortex, p is the fluid 
density, k = is the quantum of circulation, Dq{T) is 
the temperature-dependent longitudinal damping coeffi- 
cient, and Ffluc(i) is a fluctuating force with the high-T 
Markovian correlator 



(2) 



in which the angled brackets denote an average over an 
ensemble of identically prepared systems. Finally, My is 
the geometry-dependent vortex mass. For a circular trap 
geometry with a radius Uo, the vortex core radius, 

the vortex mass is 23, 24i. 
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where Ps is the condensate density and 7b is Euler's con- 
stant. 

We assume that the vortex is in a harmonic well, so 
that 



F^p{t) ^ MytvlRit) = fcoR(t) 



(4) 



in which loq is the trap frequency, and ko the 'spring 
constant'. In this case the equation of motion can 
be simplified by defining the complex position variable 



+ iR^ , the complex 'dissipation' F 



Do 



a normalized noise constant a = y^/M^,, and a nor- 
malized complex fluctuation variable £,{t), where u^(t) = 
[^fluc (*) +*-^fluc (0] ■ "The equation of motion then becomes 



— + r--M..^M..^(i). 



(5) 



This equation is most easily solved by considering the 
4-dimensional phase space position vector Q 



(6) 



In terms of Q the solution to equation ([S|) can be written. 



Q{t) = e-^*Q(0) + ae-^* / As e^' ^^^^^ 



(7) 



where Ai is the matrix 
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-1 



(8) 



From the solution [7] we may immediately calculate the 
correlator. 



{Q{t')Q\t)) =e'-^* (O(O)Qt(o)) 



(9) 



(10) 

Here a superscript dagger denotes the matrix hcrmitian 
conjugate and it is assumed that t' > t. 



ESCAPE FROM TRAP 

To get an estimate of the time it takes the vortex to 
escape from the trap we look at (|]R(t)p) as a function 
of time. Since the eigenvalues of Ai are ^(r ± A) where 
A = ■\/T^ + 4a^ is the discriminant of the characteristic 
polynomial, we see that the leading order term is for large 
t is ^ exp (2ReAt) where ReA denotes the real part of 
A. The full leading order term for (|R(t)p) in the long 
time limit is 
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2a2|A + rP 



|A|4(2i^+Re(A) 



(11) 



Note that any contributions from the initial velocity or 
position of the vortex are subdominant at large times 
t»|A|. 

Let us define the 'escape time' Tg by setting 
(|i?(re)p) = Rl, ie., the time it takes for the vortex to 
move to the edge of the circular container; we then find 



Aln 

T 



c 

K 



In 



y/K^TTplcLoRl 



1 



where we have defined 
A 

C 



ReA 

2a^|A + rp 

|A|4(2^+ReA 



(13) 
(14) 



The approximations we have made only give the leading 
order dependence of Te{Ro) on Ro, for large Ro- Three 
remarks are in order: 

(i) In two dimensions (lol| x ^ and Dq ^ T"^, so at 
low temperature we have 



(15) 
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where i? is a constant independent of Ro and Ao = 
A{T = 0) which is proportional to ^ln(i?o) for large 
Ro- 

(ii) We can also drop the inertial term completely (ie., 
let My — >• 0); this is equivalent to solving the HVI equa- 
tions ((U with an added noise term, but no inertial term. 
The exact solution for the displacement R(<) with the 
initial conditions R(0) = is then 
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which leads to the escape time 
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where C is a constant. Thus the presence of an inertial 
term in the equation of motion leads to a completely 
different temperature dependence in the escape time, 
(iii) We can compare these results to the work of Duine 

et al. [2l|, who have Te ~ In [^^^^ where G is inde- 
pendent of both T and Ro , and to that of Fedichev and 
Shlyapnikov [3], who find Te ^ T^^ In ^ j^" ^ , where 

^min the initial radial position of the vortex. We see 
that our result for t^. has a different dependence on both 
the trap radius Ro and the temperature T. The differ- 
ence in the dependence on i?o arises because the inertial 
mass of the vortex depends on the trap radius. The dif- 
ference in the temperature dependence here, compared 
to that of Fedichev and Shlyapnikov, comes (a) because 
they have no inertial mass term (b) because our stochas- 
tic force (which is absent from their equation) depends 
strongly on temperature, and (c) because Fedichev and 
Shlyapnikov have a different temperature dependence in 
Do{T) from that used here. While Duine et aZ.[2l| do 
have a stochastic force term in their equation of motion, 
they have no inertial mass term. 

Finally, we note that the most obvious way in which 
calculations of this kind can be tested is in experiments 
on 'pancake' quasi-2d Bose condensates, of the kind in- 
vestigated in MIT and Paris 14, 1^. In this context, we 
emphasize that there are several limitations to our calcu- 
lation. First, we have ignored the effect of the boundary 
of the condensate. Fetter and SvidzinskyQ introduced a 
background superflow to the equation of motion to cancel 
the normal component of the superfluid at the boundary 
of the condensate, and this will presumably affect our 
result. Second, the Thompson-Stamp derivation [lo| as- 
sumed that local deviations from the mean superfluid 
density were small. This assumption does not always 
hold for Bose-Einstein condensates, and it will be inter- 
esting to see how the results may be modified for BECs 
in the extreme compressible limit. Finally, we note that 
at sufficiently low temperatures, we expect serious de- 
partures from the modified HVI equations, coming from 



the non-local terms found by Thompson and Stamp - ex- 
perimental probes of this regime will be of great interest, 
particularly in view of the long-standing controversy over 
the correct equations of motion for a quantum vortex[3- 



CONCLUSIONS 

Using the modified form for the HVI equations that 
was found by Thompson and Stamp in the low-frequency 
semiclassical limit [lO[, we have studied the motion of a 
superfiuid vortex in a harmonic trap. We calculated the 
life-time of the vortex in a trap for this case, and find that 
the temperature dependence of the trap escape time is 
different than that found in previous calculations [s, 21 1, 
even when the inertia of the vortex is neglected. The 
vortex inertia also has a significant infiuence on this es- 
cape time. Experimental tests of results like this will be 
possible on trapped cold BECs, and will allow significant 
tests of the theory of vortex dynamics (which have been 
very difficult until now). 

This work was supported by NSERC, by CIFAR, and 
by PITP. 
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